Effective one-body dynamics in multiple-quantum NMR experiments 
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A suitable NMR experiment in a one-dimensional dipolar coupled spin system allows one to 
reduce the natural many-body dynamics into efltective one-body dynamics. We verify this in a 
polycrystalline sample of hydroxyapatite (HAp) by monitoring the excitation of NMR many-body 
superposition states: the multiple-quantum coherences. The observed effective one-dimensionality 
of HAp rehes on the quasi one-dimensional structure of the dipolar coupled network that, as we show 
here, is dynamically enhanced by the quantum Zeno effect. Decoherence is also probed through a 
Loschmidt echo experiment, where the time reversal is implemented on the double-quantum Hamil- 
tonian, Hdq oc it ^'j + ■ We contrast the decoherence of adamantane, a standard three- 

dimensional system, with that of HAp. While the first shows an abrupt Fermi-type decay, HAp 
presents a smooth exponential law. 

PACS numbers: 03.67.Pp, 03.65.Xp, 76.60.Lz, 76.90.-fd 



I. INTRODUCTION 

The new developments in nanodevices P, H, spintron- 
ics and quantum information processing critically 
rely on the control of quantum dynamics. This control is 
challenging because the manipulation of quantum states 
is crucially limited by decoherence [1, 0] • In this sense, 
much can be learned from nuclear magnetic resonance 
[1, , which offers the opportunity to tailor the interac- 
tions, and thus the time scales, and to quantify decoher- 
ence by implementing Loschmidt echoes [10]. 

The control of interaction anisotropy, e.g., the switch 
from a dipolar to an XY (planar) interaction, provides a 
tool for enhancing the transfer of quantum information 
[ill [T^ . In particular, the interactions can be sequen- 
tially turned on and off to prune some branches in real 
space so that an excitation is directed to a desired target 
through a specific pathway ,13|. By exploiting the map- 
ping between spins and fermions, spin state transfer in 
linear spin chains and rings coupled by XY interaction 
was proposed [ll] and observed in liquid-state NMR [l^ . 
Moreover, new sug gest ions that improve state transfer 
have been reported [16l [TtI. [Tst. The structurally quasi- 
one-dimensional spin systems of hydroxyapatite (HAp) 
and fluorapatite have been proposed as candidates for 
implementin g qu antum information processing in solid- 
state NMR [1^. In these systems, universal control 
has been achieved by implementing collective control to- 
gether with suitable spin manipulation at the chain ends 
|20l [2l| . All these fine control attempts might be frus- 
trated by decoherence [22, 23, 24]. Thus, the dependence 
of decoherence on nuclear spin network topology becomes 
an important issue. 
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In this work, we tailor the interactions in a one- 
dimensional (Id) dipolar coupled spin system to trans- 
form its natural many-body dynamics into effective one- 
body dynamics. The difference in dynamics is observed 
through the excitation of NMR many-body superposi- 
tion states: the multiple-quantum coherences [25|. Each 
M-quantum coherence (M-QC) collects all the superpo- 
sitions between two Zeeman states whose difference in 
total magnetic moment is the integer M. A/-QC intensi- 
ties are tested in solid-state NMR through phase codifica- 
tion techniques that allow one to follow the superposition 
weights as they are being created [2§] . 

In a homogeneous one-dimensional chain of nuclear 
spin 1/2, all spin sites have the same energy and cou- 
plings. If the spins are coupled under double-quantum 
interactions, TiDQ oc I^^j' + restricted to nearest- 

neighbors (NNs), analytical methods give closed expres- 
sions for the intensities of the multiple-quantum coher- 
ences [131 . Although Hdq acting on a thermal equilib- 
rium state excites all even-order coherences, it can be 
proved that in a one-dimensional system only zero- and 
second-order coherences are allowed [131 . The results of 
this model are compared with numerical calculations that 
include more realistic interactions and with NMR exper- 
iments in a polycrystalline sample of HAp. HAp behaves 
as a quasi-one-dimensional spin chain due to its dipolar 
coupled network structure [H, [!§] . We show that this 
anisotropy is further enhanced by a dynamical quantum 
Zeno effect (QZE). 

Decoherence is tested experimentally in HAp through a 
Loschmidt echo variant [10| based on Hdq and its rever- 
sal. The same experiment is performed in adamantane, 
a typical three-dimensional (3d) system, allowing us to 
contrast the effect of the coupling network. 

This paper is organized as follows. Section II discusses 
the multiple-quantum coherences as well as the double- 
quantum Hamiltonian. Here, the theoretical basis that 
allows one to obtain the effective one-body dynamics is 
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summarized. Section III describes the crystallographic 
and dynamical properties of HAp which make it an ef- 
fective one-dimensional system. Section IV describes the 
experimental methods. Sections V and VI, respectively, 
present numerical and experimental results for the M-QC 
dynamics. Section VII is devoted to the conclusions. 



II. MULTIPLE-QUANTUM COHERENCE AND 
EFFECTIVE ONE-BODY DYNAMICS 

In a typical solid-state NMR experiment on a system 
of N identical spins 1/2, the main interaction can be de- 
scribed by a dipolar Hamiltonian truncated with respect 
to the dominant Zeeman interaction [301 : 
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where dij = (7 ?i /(2rfj))(3cos (%) — 1) are the dipolar 
couplings, with 9ij as the angle between the internuclear 
vector Tij and the external magnetic field, and 7 as the 
gyromagnetic ratio. are the z components of the spin 
operators defined by the direction of the static magnetic 
field, and J^^and /~ are the raising and lowering opera- 
tors. In dipolar coupled spin systems at high magnetic 
field, the off-diagonal elements of the density matrix in 
the z basis, i.e., the coherences prs = can be 

labeled by the difference in the total magnetic quantum 
numbers between the states involved in the transition, 
M — rrir — nis, where \s) = nis \s), with = J^i^i- 
All the elements of the density matrix that connect two 
states whose difference in total magnetic moment is M 
contribute to the intensity of an M-QC [l^]. Although 
only single-quantum coherences {M = ±1) are directly 
observed by NMR, phase codification techniques [2^ al- 
low one to obtain information on the multiple-quantum 
coherences. 
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FIG. 1: Pathways to generate multiple-quantum coherences 
from the initial state p{0). Experimentally one follows the 
dashed arrow. This is equivalent, in a Id system with nearest- 
neighbor interactions, to the mathematical pathway indicated 
in solid arrows. 



In order to create coherences from an initial thermal 
equilibrium state, a Hamiltonian which does not com- 
mute with its density matrix is necessary. Both, the 



dipolar Hamiltonian rotated to the x axis {Tixx) and 
the double-quantum Hamiltonian {'Hdq) fulfill this re- 
quirement and are experimentally achievable: 
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In the special case of NN interactions, TLdq is unitary 
similar to the XY Hamiltonian, TixY oc IJ + Ij^ ■ 
Consequently, TLdq can simulate the TixY dynamics af- 
ter the corresponding transformation of the initial state. 
Although this relation between T-Ldq and Hxy is valid in 
one, two and three-dimensions [13, [3l|, we focus on one- 
dimensional systems, for which closed analytical results 
are available. Here, we summarize the successive trans- 
formations, developed by Doronin et al. [l^l, that enable 
this mapping. First, one applies the unitary transforma- 
tion 



U — exp(— iTr/f ) exp(— iTr/f )... exp(— ztt/j 
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to Hdq- This is a composition of tt pulses which rotate 
even- numbered spins 180° about the x axis. As a result, 
the transformed Hamiltonian is 



Hxy = UUdqU^ = (^.^ V+i + it 
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The same transformation must be applied to the initial 
state. For the thermal equilibrium state, in the high-field 
and high-temperature limit, we only consider the main 
deviation of the density matrix from the identity, which 
is the experimentally observable part, i.e., p (0) — "^^li ■ 
This leads to 



(9) 



Then, as shown schematically in Fig. [1] the dynamics 
of an initial state p(0) under H-dq is reduced to the dy- 
namics of p (0) under Hxy which, in turn, maps to a 
non-interacting fermion system [s^, [13, [IJ ■ The dynam- 
ics of this fermionic system has a closed analytical solu- 
tion when the interaction is homogeneous, c?i,i+i = d, Vi. 
Transforming back to the double-quantum dynamics, a 
closed expression for the density matrix p it) can be ob- 
tained. The intensities Jm of the M-QC are calculated 



where 



JA/(t) = Tr{pA/(t) P-M{t)}, 



PM (t) = E Prs (t) ., 
r.s 
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where restricts the sum to — rris = M. Thus, Pm 
cohects all the contributions to p due to coherences of 
order M, and p{t) = YImPm {t)- Then, the Jnit), in 
the normalized form X^a/ — 1; result in 



Jo {t) = ^Ecos^ (^Adt/hcos (^j^ 
^-W = ^Esin^(4dtAcos(^ 



(12) 
(13) 



with n ~ 1, N. This shows that only Z-QC and 2-QC 
are allowed. All other even-orders can not be created. 
Even though a closed analytical solution is not possible 
in a NN inhomogeneous case, it was shown that only zero- 
and second-order coherences are excited jllj, as what 
occurs in the homogeneous chain. 

Finally, the evolution of particular initial conditions 
po| under Id nearest-neighbor double-quantum inter- 
actions reduces to that of non- interacting ("one-ftody") 
spinless fermions. This one-body dynamics manifests 
through the presence of only 2 orders of coherence (Z-QC 
and 2-QC). We test this in Sees. V and VI by perform- 
ing numerical simulations and multiple-quantum NMR 
experiments of the dynamics under TioQ- This is con- 
trasted with the irreducible many-body dipolar dynamics 
under Hxx- 



the average over solid angle, we obtain 

^/(f^ = a(«,*))(^)'»i.5x2o. ,1.) 

where / (9, (p) is the angular function that takes into ac- 
count the angular dependence of the dipolar interaction 
and the relative orientation of the internuclear vectors 
with respect to the external magnetic field. 

There is a dynamical effect that further enhances the 
difference between these two couplings. The characteris- 
tic time for a flip-flop process within the chain is clearly 



(16) 



However, the characteristic rate of a flip-flop due to the 
weak cross-chain couplings should be estimated invoking 
the Fermi golden rule that yields [36| 



1/ ~ 1 
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and not d^/H as one might first guess. This is because 
the strong in-chain dynamics leads to an uncertainty of 
the final state over a wide excitation spectrum. Then, 
we have 



III. DYNAMICAL ENHANCEMENT OF THE 
ONE-DIMENSIONALITY BY THE QUANTUM 
ZENO EFFECT 

We perform NMR experiments in a physical system 
that behaves as a one-dimensional spin 1/2 chain. The 
system is a polycrystalline sample of hexagonal hydrox- 
yapatite, Cas (P04)3 OH, with space group P63/m. Due 
to the difference in resonance frequencies of the various 
spin nuclei, the experimental setup allows taking account 
of only the spin degrees of freedom of the ^H nuclei. The 
hydrogen spins of this sample are ordered as linear chains 
in the c direction of a hexagonal arrangement (a = b, c) 
[29| . A central chain is surrounded by six neighboring 
chains at a distance of — 9.42 A, [r^ — a). The closest 
distance between protons within a chain is — 3.44 A, 
(^n = c/2). In solid-state NMR the strongest interac- 
tion is the dipolar one. Because of the dependence of 
the dipolar couplings on the spin distance, the ratio be- 
tween the in-chain, din, and the cross-chain, dx, dipolar 
couplings for the orientation that maximizes the in-chain 
coupling is 



din 

dx 



2 X 20. 



(14) 



As we work with a polycrystal, we calculate for each chain 
orientation the ratio of the local second moment due to 
in-chain interactions, M2,in, to the local second moment 
due to the six neighboring chains, M2,x- Then, by taking 
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Equation ([TT]) states that fast in-chain dynamics makes 
already slow cross-chain dynamics even slower. This is 
a form of the QZE, which states that quantum dynam- 
ics is slowed down by a frequent measurement process 
[3?! . Spin-diffusion experiments in low-dimensional crys- 
tals showed an unexpected dimensional cross-over as a 
function of a structural parameter [38| . This cross-over 
was described as a QZE where the internal degrees of 
freedom act as measurement apparatus [1^. The con- 
cept that the measurement is played by an interaction 
with another quantum object, or simply another degree 
of freedom of the subsystem investigated, was indepen- 
dently and fully formalized b y r ecasting it in terms of 
an adiabatic theorem in Ref. [40|. It can even lead to a 
freeze of the spin swap dynamics as observed in a cross- 
polarization experiment [ll| . In the present context, Eq. 
p8| reinforces the effective one-dimensional behavior of 
HAp. 



IV. NMR EXPERIMENTAL SETUP 

The experiments were performed using a Bruker 
Avance II spectrometer operating dui a.^H resonance fre- 
quency of 300.13 MHz. We used a cross polarization 
magic angle spinning probe working in static conditions 
at room temperature with a 4mm outer diameter rotor. 
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The characterization of the dynamics of the multiple- 
quantum coherences was performed using the pulse se- 
quences shown in Fig. [21 The different orders of coher- 
ence excited under TioQ were generated using the two- 
pulse sequence shown in Fig. ^a,) [4l|, With this 
sequence, TioQ is built after a minimum number of scans 
Ns, with Ns = 2Mdos, where Mdcs is the order of coher- 
ence one desires to detect indirectly. Thus, in order to 
measure 2-QC, a minimum of four scans must be added. 
To get a better signal-to-noise ratio, the total number of 
scans must be a multiple of Ng . Therefore, the evolution 
of Mdos-QC under Hdq is built after Ng scans by adding 
signals with different phases (p. In particular, one uses 
= 0, 7r/2, TT and 37r/2 for filtering the 2-QC and (p — 0, 
7r/4, 7r/2, 3n/4, tt, 5tt/4, 3tt/2 and Tn/A for filtering the 
4-QC. In both cases, the phase of the reading pulse was 
alternated between and tt to keep only the orders of 
coherences Mdos ± nNs, with n = 0,l,2,... [H,!!!. 
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FIG. 2: (a) Selective two-pulse sequence to generate M-QC 
under an average double-quantum Hamiltonian Hdq based 
on rotations of Hzz (free evolution). It detects 2-QC (2-|-4n), 
or 4-QC (4 + 8n) intensities by applying the appropriate phase 
cyclings (see text). All pulses are of 7r/2. (b) Sequence to gen- 
erate M-QC under a rotated dipolar Hamiltonian Tixx- (c) 
Pulse sequence to generate Af-QC under a Hdq based on m 
repetitions the eight 7r/2- pulse pattern displayed in parenthe- 
sis. In (b) and (c) an free induction decay was recorded for 
each value of t and <j). The highest coherence order detected, 
rimax = 8, is governed by a digital phase shift increment, with 

A(j) = n/Umax- 

In order to encode A/-QC orders during the evolution 
under Hxx [Eq- O], we used the sequence shown in Fig. 
[Sfb) , which is a modified version of the sequence reported 
in Ref. [31 . Here, the highest coherence order detected, 
n-max = 8, is governed by the phase shift increment = 



In these sequences, the recorded free induction decays 
were the sum of 64 scans. The recycling time, Di = 3 s, 
was chosen to be longer than five times the spin-lattice 
relaxation time Ti « 500 ms. The 7r/2 pulse length was 
2.74 lis. The preparation times t, i.e., the periods evolv- 
ing under the desired effective Hamiltonian, were varied 
from 1 to 200 fis. The free evolution time ti = 0.5 /xs 
was negligible. After the mixing time and before the 7r/2 
reading pulse, a delay D = 2 ms was used to allow the 
transverse magnetization to decay. The detected signal 
was normalized to a reference FID obtained by the ap- 
plication of a 7r/2 pulse with the same number of scans. 

The two-pulse sequence used to generate Hdq was 
chosen because the fast growth of the 2-QC intensity is 
not captured with the eight-pulse sequence shown in Fig. 
djc) [23 . The last only captures a few data points in the 
time range of interest because of the minimum time of 
^ 60 /is required to accommodate the eight pulses of the 
basic unit. However, the eight-pulse sequence was ap- 
plied to implement a "Loschmidt echo" experiment, that 
is, to generate Hdq and then —Hdq, by using (j) — 0. 
We use this echo to give a measure of decoherence rates. 
In order to compare this decoherence rate in HAp with 
a widely known system, we performed the echo exper- 
iments in adamantane. Adamantane is a 3d molecu- 
lar crystal with only intermolecular dipolar interactions 
[4^ (the intramolecular interactions cancel out due to 
rapid molecular rotations). In the Loschmidt echo ex- 
periments, the preparation time was varied from 60 to 
1400 fis and the 7r/2 pulse length was 2.34 fis for HAp 
and 2.20 fis for adamantane. 

The experiments were carried out in a polycrystalline 
sample of hydroxyapatite synthesized by a modification 
of the biomimetic method reported by Zhang et al. [i^ , 
while a commercial polycrystalline sample of adamantane 
was used as provided. 



V. NUMERICAL RESULTS: 
MULTIPLE-QUANTUM DYNAMICS 

The M-QC intensities were numerically simulated us- 
ing an ensemble average of the evolution of each Zee- 
man state. The total magnetization was calculated as a 
function of preparation time t and as a function of the 
Af-QC codification phase (j>. This was obtained by evolv- 
ing each initial state under H during t and then under 
—H^, where H^ = exp(— 10/^)7^ exp(i(/>/^). Finally, a 
fast Fourier transform on (j) was applied to the magneti- 
zation to obtain the A/-QC intensities Jm (t) f25|- 

An alternative method to obtain Jm (t), which makes 
use of Eq. PH)) . was used. In this case, the pM are ob- 
tained from the elements of the density matrix calculated 
for each Zeeman state. Although this method is time 
consuming, it clearly shows how the different coherences 
contribute to the intensity of a given order. 

This second method allows us to draw some conclu- 
sions about the unitary transformations schematized in 



5 



Fig. [T] Even when TH-dq is unitary similar to TixY , an ar- 
bitrary initial condition under TIdq does not necessarily 
yield only 2 orders of coherences. In a chain with NN XY 
interaction, any excitation remains in the same subspace, 
i.e., only zero-order coherences appear. However, the 
transformed initial thermal equilibrium condition p (0) 
[Eq. ^] imposes a further restriction in the accessible 
Hilbert space in which this condition can evolve under 
XY interaction. In this case, only a portion of the ZQ- 
subspace can be reached. It is because of this restriction 
that, after transforming back to the double-quantum dy- 
namics p{t)^ only zero- and second-order coherences are 
excited. 

In order to obtain the dynamics of Jm (t) under TCdq 
and contrast this with that under Tixx , we used the first 
method described above. Since the effective Hamiltoni- 
ans 7i_DQ and TLxx are built up experimentally from the 
natural dipolar interaction, which decays with it 
becomes important to take into account the next-nearest- 
neighbor (NNN) interaction in the simulations. In a 
chain, the values of the NNN couplings are 1/8 of the NN 
ones. The simulated dynamics of the Z-QC, 2-QC and 
4-QC intensities under Jixx and 'Hdq is shown in Figs. 
13] and m for an = 10 spin chain starting at thermal 
equilibrium. Preliminary experimental results in poly- 
crystalline HAp showed that there were no detectable 
M-QC intensities after 200 /xs. Consequently, we do not 
need simulations for longer times, but we have to take 
a large enough number of spins to avoid distortions of 
the dynamics due to reflections at the chain ends. To 
verify this, we calculated the earliest time at which the 
mesoscopic echo, i.e., the revival tha t ap pears because of 
the finite nature of the system [l^, l47j |. occurs. This is 
ensured by using a single crystal at orientation 9ij = 0, 
leading to the maximum dipolar coupling dmax, which 
for HAp is dmax = 27r?i x (2.95 kHz). Any other orien- 
tation will just stretch the time scale of this curve, de- 
laying the occurrence of the mesoscopic echo. As shown 
in Fig. [21 for ten spins the mesoscopic echo appears at 
6 S/dmax ~ 325 /is. It is important to emphasize that 
by varying slightly the number of spins, the dynamics 
changes only in the neighborhood of the mesoscopic echo, 
remaining unaffected before 3.7 /i/c?max ~ 200 ^s. 

In Fig. m the 4-QC intensity dynamics in a ten spin 
chain of HAp with NNN interactions under T-Ldq and 
Tixx is displayed. In each Hamiltonian evolution, a sin- 
gle orientation of the chain and a powder average (the 
integral over solid angle of the orientation dependent dy- 
namics) were calculated. Notice that the observable we 
are using to check the effective one-body dynamics is ro- 
bust under orientation average, i.e., the non-excitability 
of the 4-QC occurs for every orientation of the chains in 
a polycrystalline sample, maintaining its null intensity. 
As can be seen, if one includes the NNN interactions in 
the chain, the intensity of the 4-QC under Hdq is not 
strictly zero. However, this intensity will not be observed 
under the typical conditions of an NMR experiment. In 
contrast, the intensity of a 4-QC under Tixx might be 
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FIG. 3: (Color online) Numerical simulations of the dynamics 
of the Z-QC, 2-QC and 4-QC intensities of HAp under Tixx 
(dashed line for 2-QC and 4-QC, and dash-dotted line for Z- 
QC) and Hdq (solid line for 2-QC and 4-QC, and dotted line 
for Z-QC) in a 10 spin chain with nnn interaction for the chain 
orientation that maximizes the coupling, dmax/?i = 27r x 2.95 
kHz. The mesoscopic echo appears at 6 ft/dmax. 



0.20 



0.00 
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FIG. 4: (Color online) Numerical simulations of the dynamics 
of the 4-QC intensities under Hxx and Hdq for a 10 spin 
chain of HAp with nnn interactions. Dashed line corresponds 
to Hxx and solid line to Hdq, both at the orientation of 
maximum dipolar coupling, dmax/fi = 27r x 2.95 kHz. The 
powder average is shown with a dotted line for < Hxx > 
and with dash-dotted line for < Hdq >■ 



observable. 

The inclusion of an extra interaction, in this case the 
NNN interaction, breaks the mapping to non-interacting 
fermions. Consequently, the system evolution is no longer 
restricted to only Z-QC and 2-QC. However, as it is 
clearly shown in Fig. [H the 4-QC under Hdq is 1 order 
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of magnitude smaller than the 4-QC excited by Tixx- 
This means that Ti-DQ still keeps the main dynamics be- 
tween Z-QC and 2-QC. Hence, one can infer that the 
effective one-body dynamics is preserved as a good ap- 
proximation. 



VI. EXPERIMENTAL RESULTS: 
MULTIPLE-QUANTUM DYNAMICS AND 
DECOHERENCE 

The pulse sequences shown in Figs. Elja) andlS^b) were 
used to generate M-QC under the effective Hamiltonians 
Ti-DQ and Tixx, respectively, from a thermal equilibrium 
state. Figure O displays the 2-QC and 4-QC intensities 
as functions of the preparation time t. There, the 4-QC 
has been multiplied by a factor of 10 because of its small 
intensity as compared with the 2-QC. While the 4-QC 
intensity under Hxx is well above the noise level, being 
evident its growth and decay, the intensity of the 4-QC 
under Hdq remains at the noise level. 



0.08 I— r 
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FIG. 5: (Color online) Dynamics of 2-QC and 4-QC intensities 
under Hdq and Ttxx in HAp implemented with the pulse 
sequence of Fig. [2] (a) and (b), respectively. Notice that the 
normalized intensities of the 4-QC are 10 times enlarged. 



In our particular Id system, the essential difference 
between the ideal Tixx and Hdq is that the first al- 
lows for the development of many orders of coherence, 
while the second allows only 2. Since higher orders of 
coherences decay at higher rates [i^ H^, we expect a 
faster decoherence in the case of Hxx- However, we 
should also assess the precision of the experimental se- 
quences used to generate these Hamiltonians. In this 
context, one should remember that our implementation 
of a multiple-quantum experiment under Hxx includes 
a dipolar Hamiltonian reversal. This involves a further 
truncation of the dipolar Hamiltonian with respect to 
the rf Zeeman interaction during the long rf pulse [30| . 



This produces additional decoherence because the trun- 
cated non-secular terms, whose magnitudes depend on 
the rf power, are not reversed [1^. The pulse sequences 
to generate Hdq may also have some limitations. It is 
known that the eight-pulse sequence produces a much 
better average Hamiltonian than the two-pulse one, espe- 
cially for long preparation times [1^ H^l ■ This is because 
the two-pulse sequence plotted in Fig. [DJa) does not av- 
erage out the chemical shift nor cancels out rf inhomo- 
geneities as the eight-pulse sequence does. For example, 
if we compare the 2-QC intensities in HAp for the eight- 
and two-pulse sequences, both of them show exponen- 
tial decay. However, the characteristic time of the first 
is Tgp « 210 /xs, while that of the second is T2p « 65 /xs, 
i.e., it is three times faster. The rapid decay of the 2-QC 
intensity with the two-pulse sequence explains the early 
occurrence of the maximum (w 15 /is) in the evolution 
of 2-QC (see Fig. [5]) as compared with the theoretical 
estimation of ft/dmax ~ 50 fis in Fig. O Because of this, 
the decay of the 2-QC for the two-pulse sequence is not 
a reliable quantifier of the decoherence of the system. 

In order to have a measure of the global decoherence 
time of the system under Hdq, we used the eight-pulse 
sequence shown in Fig. [U^c). Having minimized possible 
experimental artifacts, we expect to have a decoherence 
that reflects the properties of the sample itself (topology 
of the coupling network, defects, etc.). Following this 
idea, we compare the behaviors of HAp and adamantane 
measuring a Loschmidt echo, that is, generating Hdq 
and then —Hdq- The decays of both systems are dis- 
played in Fig. [HI The difference in the functional form 
of decay is remarkable. While a simple exponential with 
characteristic time (770 ± 50) /is holds for HAp, a 
Fermi-type curve M{t) oc 1/ [1 + exp ((< — tc)/T^)] with 
tc = (545 ± 2) IIS and — (123 ± 2) //s, provides the 
best fit for adamantane. 




FIG. 6: (Color online) Loschmidt echo experiment based on 
TLdq and —Hdq using the sequence shown in Fig[2lc) with 
= in HAp (squares) and adamantane (circles). 



It should be noticed that in adamantane, coherences 
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of very high order are generated quite rapidly. Indeed, 
coherences of orders M > 100 are well defined after 0.5 
ms [H, [5l|, indicating the huge portion of the Hilbert 
space explored through TioQ in this system. As shown 
by the Fermi-type curve, the coherence of such highly 
interacting system is not sustained beyond a critical time 
tc where a sort of "catastrophe" seems to occur. A similar 
behavior is observed in simulations of highly interacting 
systems, either fermions or bosons, whose coherence also 
decays following a Fermi- type curve [H, HI]. In those 
works, a self-consistent approximation allows one to see 
this critical stage as the triggering of a nonlinear loop. 

In contrast with adamantane, the decoherence of HAp, 
as seen from the Loschmidt echo, occurs smoothly fol- 
lowing an exponential law. This sort of decay has been 
seen in chaotic one-body systems in semiclassical states 
where the perturbation effects are limited [s^, [H, [s^ . 
Hence, this decay is consistent with the restricted dy- 
namics imposed by the low connectivity of a Id system. 
Furthermore, as the dominant dynamics is that of the 
non-interacting fermions, the residual interactions and 
the experimental imperfections define the "environment" 
that produces the exponential decoherence. 

Although the observed decoherence rate seems to be 
somewhat fast to enable a straightforward quantum in- 
formation application, the exponential decay of the dy- 
namics of the Id system may be easier to manipulate 
than the dynamics of the 3d system. On the other hand, 
the 3d system presents a short-time behavior that could 
be nicely exploited to implement quantum operations, 
because the coherence is lost at a very low rate. Fur- 
ther experimental designs are necessary to confirm the 
origin of these different functional forms, and to quan- 
tify the factors determining the respective characteristic 
decay times in HAp and tc and r^p in adamantane. 

VII. CONCLUSION 

We have shown that the Af-QC intensities under a 
double-quantum Hamiltonian in HAp behave as effective 
one-body dynamics. This has been observed through sev- 



eral experiments where the evolutions of the intensities 
of the 2-QC and 4-QC were studied under the action 
of Hdq- These results were contrasted with the many- 
body dynamics induced by TLxx- No coherence orders 
above 2 appear under H-dq, while they do appear un- 
der TLxx- In both cases, the dynamics remains mainly 
one-dimensional as the natural anisotropy of HAp is en- 
hanced by the quantum Zeno effect. 

The global decoherence of HAp under Hdq was com- 
pared with that in adamantane, a regular 3d system, 
whose genuine many-body dynamics is manifested by the 
rapid excitation of very high orders of coherence. The 
coherence decays in both systems follow completely dif- 
ferent functional forms. 

In summary, we have addressed two main points: 

(1) We confirmed the mapping of a nearest-neighbor 
one-dimensional spin system under a double-quantum in- 
teraction to a non interacting fermion system. This map- 
ping was tested through one of its main consequences: 
the non excitability of 4-QC under Hdq- 

(2) We evaluated the decoherence through a Loschmidt 
echo experiment based on a double-quantum Hamilto- 
nian. The restricted dynamics induced by the low con- 
nectivity space leads to the appearance of a smooth ex- 
ponential decoherence, while the dynamics in a high con- 
nectivity space shows a sudden drop in coherence. 

These results indicate that, in spite of residual inter- 
actions, HAp can be used as a "quantum simulator" for 
non interacting fermion dynamics. 



Acknowledgments 

The authors thanks R. H. Acosta, Y. Garro Linck and 
A. K. Chattah for fruitful experimental discussions, as 
well as L. E. F. Foa Torres and F. Pastawski for sug- 
gestions. This work was made possible through the fi- 
nancial support from CONICET, ANPCyT and SeCyT- 
UNC. PRL and HMP acknowledge the hospitality of the 
Institute for Materials Science TU-Dresden and the MPI- 
PKS. 



[1] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 
A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, Science 309, 2180 (2005). 

[2] M. V. G. Dutt, L. Childress, J. Liang, E. Togan, J. Maze, 
F. Jelezko, A. S. Zibrov, P. R. Hemmer, and M. D. Lukin, 
Science 316, 1312 (2007). 

[3] D. D. Awschalom and M. E. Platte, Nat. Phys. 3, 153 
(2007). 

[4] D. P. DiVincenzo, Science 270, 255 (1995). 

[5] K. C. Nowack, F. H. L. Koppens, Y. V. Nazarov, and 

L. M. K. Vandersypen, Science 318, 1430 (2007). 
[6] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 
[7] W. H. Zurek, F. M. Cucchietti, and J. P. Paz, Acta Phys. 



Pol. B 38, 1685 (2007). 
[8] D. G. Cory, A. F. Fahmy, and T. F. Havel, Proc. Natl. 

Acad. Sci. U.S.A. 94, 1634 (1997). 
[9] L. M. K. Vandersypen and I. L. Chuang, Rev. Mod. Phys. 

76, 1037 (2005). 
[10] R. A. Jalabert and H. M. Pastawski, Phys. Rev. Lett. 

86, 2490 (2001). 
[11] G. A. Alvarez, E. P. DanieU, P. R. Levstein, and H. M. 

Pastawski, J. Chem. Phys. 124, 194507 (2006). 
[12] G. A. Alvarez, E. P. DanieU, P. R. Levstein, and H. M. 

Pastawski, Phys. Rev. Lett. 101, 120503 (2008). 
[13] G. A. Alvarez, PhD thesis, Universidad Nacional de 

Cordoba (unpubUshed) e-print arXiv:0705.2350 (2007). 



8 



[14] H. M. Pastawski, G. Usaj, and P. R. Levstein, Chem. 

Phys. Lett. 261, 329 (1996). 
[15] Z. L. Madi, B. Bmtshor, T. Schultc-Hcrbriiggon, 

R. Briischweiler, and R. R. Ernst, Chem. Phys. Lett. 268, 

300 (1997). 

[16] S. Rose, Phys. Rev. Lett. 91, 207901 (2003). 

[17] M. Christandl, N. Datta, A. Ekert, and A. J. Landahl, 

Phys. Rev. Lett. 92, 187902 (2004). 
[18] A. Kay, Phys. Rev. Lett. 98, 010501 (2007). 
[19] T. D. Ladd, J. R. Goldman, A. Dana, F. Yamaguchi, and 

Y. Yamamoto, arXiv:quant-ph/0009122v4 (2001). 
[20] P. Cappellaro, C. Ramanathan, and D. G. Cory, Phys. 

Rev. Lett. 99, 250506 (2007). 
[21] P. Cappellaro, C. Ramanathan, and D. G. Cory, Phys. 

Rev. A 76, 032317 (2007). 
[22] H. M. Pastawski, P. R. Levstein, G. Usaj, J. Raya, and 

J. Hirschinger, Physica A 283, 166 (2000). 
[23] P. R. Levstein, G. Usaj, and H. M. Pastawski, J. Chem. 

Phys. 108, 2718 (1998). 
[24] H. G. Krojanski and D. Suter, Phys. Rev. A 74, 062319 

(2006). 

[25] J. Baum, M. Munowitz, A. N. Garroway, and A. Pines, 
J. Chem. Phys. 83, 2015 (1985). 

[26] R. R. Ernst, G. Bodenhausen, and A. Wokaun, Prin- 
ciples of Nuclear Magnetic Resonance in One and Two 
Dimensions (Oxford University Press, Oxford, 1987). 

[27] S. L Doronin, I. L Maksimov, and E. B. Fel'dman, Zh. 
Eksp. Teor. Fiz. (2000). 

[28] G. Cho and J. P. Yesinowski, Chem. Phys. Lett. 205, 1 
(1993). 

[29] G. Cho and J. P. Yesinowski, J. Phys. Chem. 100, 15716 

(1996). 

[30] C. P. Slichtcr, Principles of Magnetic Resonance 

(Springer- Verlag, New York, 1992). 
[31] C. A. Perez-Delgado, M. Mosca, P. Cappellaro, and D. G. 

Cory, Phys. Rev. Lett. 97, 100501 (2006). 
[32] E. H. Lieb, T. Schultz, and D. C. Mattis, Ann. Phys. 

(N.Y.) 16, 407 (1961). 
[33] E. B. Fel'dman and M. Rudavets, Chem. Phys. Lett. 311, 

453 (1999). 

[34] E. P. Danieh, H. M. Pastawski, and P. R. Levstein, Chem. 

Phys. Lett. 384, 306 (2004). 
[35] S. L Doronin and E. B. Fel'dman, SoUd State Nucl. 

Magn. Reson. 28, 111 (2005). 



[36] E. Rufoil-Fiori and H. M. Pastawski, Chem. Phys. Lett. 

420, 35 (2006). 
[37] B. Misra and E. C. G. Sudarshan, J. Math. Phys. 18, 

756 (1977). 

[38] P. R. Levstein, H. M. Pastawski, and R. Calvo, J. Phys. 

Condens. Matt. 3, 1877 (1991). 
[39] H. M. Pastawski and G. Usaj, Phys. Rev. B 57, 5017 

(1998). 

[40] P. Facchi and S. Pascazio, Phys. Rev. Lett. 89, 080401 
(2002). 

[41] M. A. Voda, D. E. Demco, J. Perlo, R. A. Orza, and 
B. Bliimich, J. Magn. Reson. 172, 98 (2005). 

[42] R. H. Acosta, G. A. Monti, D. A. Vega, M. A. Villar, and 
E. M. Valles, Macromolecules 39, 4788 (2006). 

[43] M. Munowitz and A. Pines, Adv. Chem. Phys. 66, 1 
(1986). 

[44] H. Cho, T. D. Ladd, J. Baugh, D. G. Cory, and C. Ra- 
manathan, Phys. Rev. B 72, 054427 (2005). 
[45] J. Baum and A. Pines, J. Am. Chem. Soc. 108, 7447 

(1986). 

[46] Q. Zhang, Y. Leng, and R. Xin, Biomatcrials 26, 2857 
(2005). 

[47] H. M. Pastawski, P. R. Levstein, and G. Usaj, Phys. Rev. 

Lett. 75, 4310 (1995). 
[48] H. G. Krojanski and D. Suter, Phys. Rev. Lett. 93, 

090501 (2004). 

[49] C. M. Sanchez, H. M. Pastawski, and P. R. Levstein, 

Physica B 398, 472 (2007). 
[50] M. Munowitz, Coherence and NMR (John Wiley & Sons, 

New York, 1988). 
[51] C. M. Sanchez, P. R. Levstein, R. H. Acosta, and A. K. 

Chattah, unpublished (2008). 
[52] G. Manfredi and P.-A. Hcrvioux, Phys. Rev. Lett. 97, 

190404 (2006). 

[53] G. Manfredi and P.-A. Hervieux, Phys. Rev. Lett. 100, 
050405 (2008). 

[54] F. M. Cucchietti, H. M. Pastawski, and R. A. Jalabert, 

Phys. Rev. B 70, 035311 (2004). 
[55] T. Gorin, T. Prosen, T. Sehgman, and M. Znidaric, Phys. 

Rep. 435, 33 (2006). 
[56] P. Jacquod and C. Petitjean, Adv. Phys. (to be pub- 

hshed) e-print arXiv:0806.0987vl (2008). 



